Mott-Hubbard Scenario for the Metal-Insulator Transition in the Two Dimensional 

Electron Gas 



O 
O 

(N 

< 



■4— » 
•4-* 



c 

O 

o 



> 
o 

o 

oo 
O 

: 

o . 
-i— > ■ 

s; 

T3 ■ 

: 

o . 
o ■ 



X 



Ping Sun 1 

1 Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854-8019 

(Dated: February 6, 2008) 

By comparing the responses to an in-plane magnetic field near the metal-insulator transition 
(MIT), we find that the observed MIT in Si MOSFETs can be described by the non-perturbative 
Mott-Hubbard scenario. Interrelations between independent measurables are uncovered and con- 
firmed by replotting the experimental data. A universal critical energy scale vanishing at the MIT 
is extracted from the experimental data and the critical exponent found. 
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The properties of an electron gas (EG) are controlled 
by the average inter-electron distance. As the elec- 
tron density is tuned from high to low, the Coulomb 
interaction, favoring localization of the electrons, be- 
comes dominant over the kinetic energy and the system 
evolves from a metal to an insulator |l|. Puzzling mag- 
netotransport experiments on the silicon metal-oxidc- 
semiconductor field-effect transistors (Si MOSFETs) in 
the past ten years00] have reactivated broad interest in 
the subject because of possible relations to the detection 
of a quantum critical point (QCP), the realization of the 
Wigner crystal[l|, and the interplay between correlation 
and disorder ^ |j, 0, Q. By comparing the responses to 
an in-plane magnetic field near the metal-insulator tran- 
sition (MIT), we find that the observed MIT is controlled 
by the non-perturbative Mott-Hubbard (M-H) critical 
end point at a low but non-zero temperature [a. 1^. llCj . 
instead of a QCP. This also contradicts the long-standing 
speculation that the non-local part of the Coulomb inter- 
action plays an essential role in the MIT, as it does in the 
Wigner crystalizationQ. The M-H scenario predicts in- 
terrelations between independent measurables which are 
confirmed by replotting the experimental data. The un- 
covered sample independent behavior indicates that the 
disorder does not play an active role. A critical energy 
scale is extracted from the experimental data and the 
critical exponent found. 

An MIT in the two dimensional (2D) EG was first ob- 
served via the conductivity measurements in ultra clean 
Si MOSFETs[ll|. When the temperature is decreased 
below ~ 2K, the conductivity is observed to increase (de- 
crease) monotonically at densities above (below) a cer- 
tain non-universal critical value n c ~ 10 11 cm -2 . In this 
temperature range, the conductivity curves in a given 
phase can be all collapsed into a single curve by rescal- 
ing the temperature|l2l|. The observation clearly points 
to an unidentified type of critical behavior. It was also 
found 01 that the scaling behavior is sample independent 
down to the lowest temperatures reached experimentally 
(< 100mA'). The most likely driving force is then the 
interaction induced correlation, instead of the disorder. 

In the subsequent experiments in Si MOSFETs, it was 
found that near n c applying an increasing in-plane mag- 



netic field Hu causes the conductivity first to decrease 
and then, above a certain field H ai saturate at a roughly 
field independent value [l3L fTi!] . As the temperature is 
lowered, H a first decreases linearly and then approaches 
a constant. This constant depends on the electron den- 
sity and vanishes as the transition is approached from 
the metallic side^|. At ifn > H a , the conductivity is 
always insulating like, decreasing as the temperature is 
lowered, irrespective of whether the density corresponds 
to a metallic or an insulating phase in the absence of the 
field^{|. Shubnikov-de Haas (SdH) measurements [l6j 
show that the oscillation period is halved when Hu > H a . 
This has been taken as an indication that the entire EG 
is fully polarized when the conductivity saturates in the 
fieldhj. Since H a vanishes at n = n c and T = 0, this 
could point to a Stoner instability^]}, although it would 
be hard to understand the almost coincidence of the 
Stoner instability with the MIT[i3|. Direct [H EU and 
indirect jFJ Hj, |2Q| measurements show that the homo- 
geneous spin susceptibility is enhanced in the approach 
to the transition and is related to a mass enhancement 
since the g-factor changes little|2l|. 

The perturbative scaling theory for non-interacting 
disordered electron systems shows that no true metal- 
lic phase exist in 2D due to the disorder Q. How- 
ever, the Coulomb interaction competes with the local- 
ization and stablizcs the metallic phased H 0- The 
relevance of the scaling theory to the experiments was 
checkedj^. Starting from the Fermi liquid, pertur- 
bations were used to study the temperature depen- 
dence of the conductivity [22l l23l | and the strong mass 
enhancement p3| . While they describe the 2DEG to cer- 
tain extent, these perturbations are not justified since the 
experiments 0,0] on Si MOSFETs show that the transi- 
tion happens when R$, the average inter electron dis- 
tance measured in units of the Bohr radius, is > 10. The 
significance of this is that Rs is also given by the ratio 
of the interaction energy, i.e. the perturbation, to the 
kinetic energy. 

The only non-perturbative, correlation-driven MIT we 
understand is the M-H transitionl8tl9 Ul0l[ . The M-H sce- 
nario was related to the 2DEG beforepM l2r| . However, 
the clear connections between the two become evident 
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only through studying the magnctotransport behavior. 
Within this scenario, we manage to explain the exper- 
iments and uncover previously unnoticed interrelations 
between independent observables. 

The M-H transition, driven by the ratio of the on-site 
interaction U to the bandwidth W, is best revealed by the 
behavior of the local density of states (LDOS) as shown 
by the dynamical mean field theory (DMFT)JlfJ. Near 
the transition and on the metallic side, the LDOS has 
a three peak structure (Fig. 1), with a central peak at 
the Fermi level and two side bands separated in energy 
by roughly U. The former consists of extended quasipar- 
ticle states; its width is proportional to the quasiparti- 
cle weight Z. The latter, the upper and lower Hubbard 
bands, consist of localized states. In between the peaks, 
the LDOS is featureless, corresponding to an incoherent 
background. As one moves towards the insulating side by 
increasing U /W, the width of the central peak decreases 
to zero. DMFT found[H| that, irrespective of the bare 
dispersion details, the MIT happens at U/{W/2) ~ 3 and 
is first-order at zero temperature. The first-order transi- 
tion line extends to non-zero temperatures and ends in a 
critical point at T c ~ 0.05 • (W/2) beyond which the tran- 
sition becomes a crossover. This picture was confirmed 
experimentally only recently, e.g. in the photoemission 
studies|13 of V 2 3 at T > T c . At T < T c the physics is 
usually non-universal and depends on the system details. 

We argue that the magnetotransport experiments on 
the 2DEG can be explained by the evolution in a mag- 
netic field of the three peak LDOS. On the metallic side, 
an increasing in-plane magnetic field first splits the quasi- 
particlc peak while reducing its height. This reduces the 
LDOS near the Fermi surface and thus the conductivity. 
When the field is strong enough so that the quasiparti- 
cle peak is either entirely split apart or suppressed, the 
LDOS left over near the Fermi surface derives only from 
the incoherent background. Consequently the conduc- 
tivity saturates at a value determined by the short mean 
free time of those states and its temperature dependence 
becomes insulating-like. The field at which this satura- 
tion appears is H a oc (1/<j/xb) - max (ZW, ksT), meaning 
the Zeeman energy should be strong enough to split or 
suppress the quasiparticle peak and overcome the ther- 
mal excitations. So as T is reduced, H a first decreases 
linearly and then saturates to a constant cx Z . Since Z re- 
duces to zero as the transition is reached from the metal- 
lic side, H a (T = 0) vanishes accordingly. Meanwhile, the 
effective mass (cx l/Z) gets enhanced and diverges. 

To support the above explanation, we have computed 
the evolution in an in-plane magnetic field of the LDOS 
near the M-H transition using a half-filled lattice Hub- 
bard model, 

H = Y [ -^L (e k + aH») C\ C r 



where we used Ni t(T = G irj Ci^ and set the (/-factor = 2. 
We approximate the first Brillouin zone by a disk. The 
momentum cut-off A = 2y / 7r ensures the normalization of 
the momentum integral and sets the lattice constant to be 
27r/A. We use the free dispersion ek = k 2 /2 specific to the 
EG. The bandwidth is thus W = 2n. This effective model 
is valid when the thermal and Zeeman energies are small 
comparing to the bandwidth. We mention, though, even 
at low temperatures there is a non-zero possibility that an 
electron jumps to a higher energy state. In terms of the 
one-band lattice model, this is like an electron hopping 
away from the lattice sites and into the interstitials. 

In this dimensionless model for the 2DEG, the en- 
ergy is measured in units of Eq = 2esi?~ 2 , where eb = 
mbe 4 /(2e 2 h 2 ), with the electron band mass mj, ~ 0.2m e 
and the dielectric constant for the Si MOSFET e ~ 7.7. 
Since in the experiments, the MIT is observed around 
R s ~ 10, we estimate Eq to be ~ 10.6AT. The cor- 
responding magnetic field unit is, using g = 2, Hq = 
Eq/ (nBTTie/mb) ~ 3. IT. A rough estimation of the non- 
universal T c for the M-H critical end point is possible at 
this stage. Applying the DMFT estimation in the current 
model, we obtain T c ~ 1.7 K. The mean field estimation 
is about one order of magnitude higher than the experi- 
mental value T c < IOOjtiA'. 

To c omp ute the conductivity, we use the Kubo 
formulajlOj. In the EG, the current operator is given 
by M ~ EkJk + q/2)Cl +? C % a . Only the particle- 
hoel bubble survives in the current-current correlation in 
DMFT [13. The conductivity in unit of e 2 /(hd), 

,. 1 /" d 2 k 1 v-^ 

k 2 G a (k,ip n )G a (k,ip n +iu), 

with p n = (2n + \)tt / f3. 

We solve this model numerically via DMFT, u sing 
quantum Monte Carlo (QMC) as the impurity solver flflj. 
In Fig. 1 the evolution of the LDOS is presented at an 
inverse temperature f3 = 5.0 and an on-site interaction 
U = 7.0. The latter is chosen so that the system is on 
the metallic side and close to the MIT. The conductivity 
saturation behavior is shown in the inset of Fig. 2. The 
T vs H a relation is shown in the main plot. In the tem- 
perature region solved, H a depends linearly on T, similar 
to that observed experimentally 01 . 

Although it provides us with some essential features 
regarding the MIT in 2DEG, this effective model con- 
tains certain limitations. These include the lack of the 
incoherent hopping of the electrons into the interstitial 
space. This should be responsible for the observed sat- 
urated conductivity being non-zero at T > 0, instead of 
vanishing as shown in the inset of Fig. 2. Actually it 
was found experimentally that the saturated conductiv- 
ity increases with the temperature and at T > 1AT it can 
reach the same order of magnitude as that measured at 
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FIG. 1: The evolution of the LDOS in an external in-plane 
magnetic field H\ \ . At zero field, the LDOS shows a quasi- 
particle peak at the Fermi energy (E = 0) together with two 
side bands. An increasing H\\ causes the peak to split and 
finally get suppressed. At Hu ~ 0.7 when the quasiparticle 
peak is fully split, the LDOS shows that the localized spins 
are still not entirely polarized. In this model calculation, the 
strong incoherent background as observed experimentally in 
the M-H systems [27I is absent from the LDOS. 



a metallic density and with no external field [T^. SdH 
measurements then become possible even at H\\ > H a . 
The observed halving of the SdH oscillation period when 
H 1 > H a is thus due to the full polarization of the itin- 
erant electrons, including those within the quasiparticle 
peak and those involved in the interstitial hopping. The 
electrons localized in the lower Hubbard band do not 
contribute to the transport. 

While only one energy scale (Z) controls the critical 
transport behavior, two distinct energy scales are re- 
sponsible for the spin polarization, the super exchange 
(oc l/U) for the local moments and the quasiparticle peak 
width (oc Z) the itinerant electrons. Near the transition, 
Z — > and it costs almost no energy to polarize the 
itinerant electrons. Meanwhile, a non-zero exchange en- 
ergy must be overcome to polarize the local moments. 
This behavior can be seen in the LDOS at H ~ 0.7 
in Fig. 1, where the quasiparticle peak is already fully 
split while the local moments are only partially polar- 
ized. A direct proof of the two-energy-scale behavior 
could be found from a measurement of the magnetiza- 
tion and a study of its saturation behavior in a magnetic 
field[2^|. For the 2DEG, no consensus has been reached 
experimentally^^ ^| , although indications favoring the 
existence of the local moment band were rer>orted|l8j|. 
This convincing experiment would allow to discern the 
M-H scenario from the Stoner picture. In the latter there 
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FIG. 2: The conductivity calculated via DMFT[T(1 as a func- 
tion of the magnetic field at different temperatures is shown 
in the inset. When the field is strong enough, the conductiv- 
ity saturates. The saturation field H a is determined when the 
conductivity reaches a = 2 x 1CP 3 ire 2 /h, which is about the 
order of the QMC statistical error. (The choice of this value 
does not change the qualitative behavior of H a .) vs T is 
shown in the main plot by the circles, together with a linear 
fitting by the line. The intersection of the line at T = gives 
an estimation ZW ~ 0.75 at U = 7.0. 
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FIG. 3: Replotting of experimental data in terms of Z vs 
Sn = n/ric — 1. Triangles correspond to the rescaling tem- 
perature To on the metallic side of sample Si- 12a reported in 
Ref . [l2^ . It is transformed as Z = T^ 2 ' 5 /10.5. The diamonds 
correspond to H a (T = 0) measured in Ref. ll3^ . It is rescaled 
as Z = H a (T — 0)/3.9. The circles represent the effective 
mass reported in Ref. |2l| which is plotted as Z — mb/m*. 
The values of n c are taken from the corresponding papers. 



Close to 5n = 0, Z can be fitted to Z oc Sn 



This relation 



is valid up to densities almost three times the critical density. 



is only one magnetic energy scale and all the electrons 
are polarized simultaneously by a vanishing field near 
the transition. 

According to the M-H scenario, all the critical behav- 
iors are related to Z. The saturation magnetic field, 
H a (T = 0) oc Z should vanish and the effective mass 
to* oc 1/Z diverge in approaching the transition from 
the metallic side. Besides, it is also responsible for the 
vanishing energy scale Tq as revealed by the universal 
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scaling of the conductivity [T^ . so To <x Z a . In Fig. 3, 
we transformed, rescaled and replotted the experimen- 
tal scaling temperature To data^^i the H a datap^|. 
and the m* data[2l]]. Although belonging to indepen- 
dent measurements on different samples, the H a and m* 
data sets collapse well when their densities overlap, re- 
vealing the previously unnoticed interrelations that we 
derived from the M-H scenario. From Fig. 3, we find 
that Z (x Sn 061 for Sn < 2. This universal relation can 
be used to calibrate the critical behaviors. Note that, 
the DMFT solution gives Z oc Sn, since[I3 Z oc U c - U 
and the mapping from the density to U is expected to 
be analytic even around the transition. The deviation of 
the critical exponent is not unexpected for a mean field 
theory. 

The M-H picture captures the universal critical be- 
havior when T c < T < 2K, where the conductivity scal- 
ing was achieved and all the other above-mentioned 
behaviors were observed 0, Q The non-universal prop- 
erties specific to the 2DEG at T <C T c were previously 
studied. Monte Carlo simulations [2^ showed, in between 
the paramagnetic metallic and the Wigner crystal phases, 
there lies at least one other phase, the polarized liquid 
phase. A similar result was obtained by studying the 
free energy. It was shown that a direct transition is for- 
bidden and various intermediate phases are present (^3]. 
At T > 2K the system crosses over to another differ- 
ent behavior. It was observed that the conductivity as 
a function of temperature is insulating like|l2j|. irrespec- 
tive of the density being above or below n c . This is likely 



due to the Anderson localization behavior [4J before the 
correlation sets in and the M-H critical point takes over 
at lower temperatures. 

We have managed to establish that the non- 
perturbative M-H critical point describe the observed 
MIT in 2DEG. Our theory provides a framework for un- 
derstanding the interrelations between the independent 
measurables, which is confirmed by replotting the exper- 
imental data. A universal energy scale vanishing at the 
transition is uncovered. In the light of the M-H scenario, 
we suggest further experiments on the 2DEG of measur- 
ing the magnetization saturation behavior in an in-plane 
magnetic field. This will allow to reveal the predicted 
two-energy-scale behavior near the transition. Subjects 
remain to be clarified including the proper descriptions of 
the incoherent hopping into the intcrstitials and the inef- 
fectiveness of the non-local part of the partially screened 
Coulomb interaction. To fulfill these purposes, a con- 
trolled derivation of an effective low energy lattice model 
for the EG will be helpful. 
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